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ABSTRACT 

We have obtained medium-resolution infrared K band spectra of 44 giants in seven heavily reddened 
clusters in the Galactic bulge, as well as 12 giants in oj Centauri. We measure the equivalent widths 
of the Na doublet, the Ca triplet, and the CO band head, and then apply the new technique of Frogel 
et al.to determine the metallicity of each star. Averaging these values, we estimate the metallicity for 
each cluster, and compare our new [Fe/H] values with previous determinations from the literature. 
Our estimates for each cluster are: NGC 6256 (-1.35), NGC 6539 (-0.79), HP 1 (-1.30), LiUer 1 
(-0.36), Palomar 6 (-0.52), Terzan 2 (-0.87), and Terzan 4 (-1.62). We briefly discuss differences 
between the various [Fe/H] scales on which it was possible to base our calibration, which is found to 
be the largest uncertainty in using this technique to determine metallicities. 

Subject headings: (Galaxy:) globular clusters: individual (NGC6256, NGC6539, HPl, Lillerl, Palo- 
mar6, Terzan2, Terzan4, omega Centauri) 



1. INTRODUCTION 

The globular clusters of the Milky Way consist of 
two distinctly separated groups, each with approxi- 
mately Gaussian metallicity ([Fe/H]) distributions. This 
bimodal ity w as first quan tified bv [ Fr eeman fc NorrisI 
(CMIl), IZhml Emh . and lArmandroff fc ZinnI il98ar 
The mean [Fe/H] values for the two groups are approx- 
imately —0.5 and —1.6. As detailed by these authors 
and many others since then, the bimodality is not con- 
fined just to the distribution over [Fe/H]. The two groups 
have quite distinct spatial distributions and kinematical 
properties as well. The metal-poor group contains clas- 
sical, population II halo objects in terms of their loca- 
tion and space motions, in addition to their low metal- 
licities. Clusters in the metal-rich group, on the other 
hand, are mostly confined to the same volume of space 
as the Galactic bul ge. Of the 33 clusters with [Fe/H] 
> -0.7 (data from iHarrij l|1996D as updated in 2003 
February), three-quarters of them lie within 20° of the 
Galactic center, and half lie within 10° of the center. Not 
surprisingly, these clusters are also heavily reddened; the 
average E{B — V) of the sample that lies within 10° 
of the Center is 1.41, which is equivalent to an average 
visual extinction {Av) of 4.4 mag. The problems asso- 
ciated with optical studies of these clusters are further 
compounded by the fact that many of them lie in very 
crowded fields and are not themselves of high surface 
brightness or concentration. 

An obvious question to ask concerning the metal-rich 
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globular clusters is whether they are in some way physi- 
cally associated with the bulge, sharing a common origin 
and evolution, as has been argued by many authors (e.g. 
Minniti 1995b; Cote 1999); or is the similarity in spatial 
distribution of the clusters and bulge stars merely coinci- 
dental and due simply to the influence of the same grav- 
itational potential. If the former is the case, then knowl- 
edge of these metal-rich clusters would be particularly 
relevant in enhancing our understanding of the formation 
and evolution of the bulge itself. Accurate abundance de- 
terminations for these clusters would be of particular im- 
portance. However, because of the difficulties associated 
with observing these clusters in the optical, spectroscopic 
abundances based on observations of individual stars are 
available for only a small fraction of those that lie closer 
than 15° from the Galactic center. The situation is some- 
what better in terms of deriving physical parameters 
from optical color-magnitude diagrams (CMDs); of par- 
ticular note is the prog ram being carried out by Ortolani 
and his collaborators iOrtolani. Bica fc Barbuvl Il997ai 
and references therein). 

In working toward this goal of learning more about 
these heavily reddened bulge clusters, we have recently 
developed a new method for the determination of [Fc/H] 
for globular clusters based on the measurement of ab- 
sorption features due to Ca, Na, and CO in the K band 
spectra of a few o f the brightest giants in each cluster 
ijFrogel et al.ll200ll hereafter Paper I). The practical ad- 
vantages of the near-IR technique we have developed are 
several and important. The necessary spectroscopic ob- 
servations are at a resolution of only « 1500 and need 
only be obtained for the brightest half-dozen or so giants 
in a cluster. Knowledge of cluster reddening or distance 
is not needed, although if dereddened colors and abso- 
lute K band magnitudes are known, the accuracy of the 
resulting [Fe/H] value is further improved. Furthermore, 
because of the close coupling between [Fe/H] and the lo- 
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cation and morphology of a cluster's red giant branch 
(RGB) , our technique is largely insensitive to the choice 
of stars so long as they are on the upper 3 or 4 mag 
{K band) of the RGB. This advantage was also demon- 
strated in our study of gian ts in the field of the Galactic 
bulge (Ramirez et alJ 1200(1) . In the K band the giants 
near the top of a cluster's RGB are very bright in con- 
trast to their relative faintness in V, which is due to 
severe molecular blanketing. And, finally, a consequence 
of the previous fact is that the effects of crowding and 
confusion with field stars are significantly reduced. 

2. SELECTION OF CLUSTERS, OBSERVATIONS, DATA 
REDUCTION 

For this paper we chose seven poorly studied globular 
clusters located in or near the Galactic bulge, which are 
believed to have [Fe/H] > — 1.0 and/or are heavily red- 
dened. For five of these clusters, our work is the first 
detailed spectroscopic study of individual cluster stars. 
Table n] summarizes some of the properties of these clus- 
ters. 

Most previous estimates of [Fe/H] values for these 
seven clusters are from the characteristics of their op- 
tical or near-IR GMDs, or from integrated spectroscopy. 
Integrated spectroscopy of globular clusters that are in 
the very crowded Galactic bulge region, especially opti- 
cal spectroscopy, will almost always tend to overestimate 
[Fe/H] because of contamination from metal- rich bulge 
giants (see the discussion in Sectional. 

Our spectroscopic data were obtained on the Blanco 4 
m telescope at the Cerro Tololo Inter- American Obser- 
vatory (GTIO) wi th either GTIO's n ear-IR spectrome- 
ter (IRS; R=1650. lDePov et al.llT990|) or with the Ohio 
State InfraRed Ini ager/Spectrograph (OSIRIS; R=1380, 
iDePov et aLlllOQ^ . Spectral coverage was generally be- 
tween 2.17 and 2.34/im. Most of our observations were 
made on the same nights as those for the calibration 
clusters (Paper I); details of the observations and data 
reduction for these seven clusters are identical to those 
described in Paper I, and so will be only briefly described 
here. 

Each spectrum begins as a set of ~ 10 individual spec- 
tra taken by stepping the star along the slit to eliminate 
the effects of uneven slit illumination, bad pixels, and 
fringing. The slit is several arc minutes long, which per- 
mits simultaneous star and sky measurements. The ar- 
ray response signature is first removed by dividing each 
of the individual frames by a dome flat field. The sky 
frame for each star is formed by median-combining the 
individual frames in each set, thus removing all traces of 
the star itself, and then subtracting it from each spectral 
frame. Each spectra in a set is then extracted using the 
IRAF APALL routine. The final spectrum for a star is 
the average of the extracted individual spectra. 

The wavelength calibration for each star is determined 
by fitting to approxi mately 12 atmospheric OH lines 
l|01iva fc OriglialflQfll visible on a sky spectrum, where 
the sky spectrum is extracted from the sky frame us- 
ing the program star as an extraction template. We fit 
the observed wavelengths to the laboratory ones with a 
second order Legendre polynomial. 

To remove telluric absorption features from each fi- 
nal stellar spectrum, we divide them by a normalized 
standard star spectrum to correct for telluric absorp- 



tion features and then multiply by a 10,000 K blackbody 
spectrum to restore the over all shape of the spectrum. 
These standard stars are of early A spectral type or hot- 
ter, with no significant spectral features in the observed 
wavelength regime, except for Brackett 7 at 2.16 /im. 
They are observed as close in time and air mass as prac- 
tical to the observations of each cluster star. 

The technique we use to measure the equivalent widths 
(EWs) of the spectral features - EW(Na), EW(Ga), and 
EW(CO) - is described in Paper I. Table[21gives the con- 
tinuum and line regions over which we integrate. Note 
that EW(GO) does not include any contribution from 
the i^CO band head. EW(Ca) and EW(Na) have con- 
tinuum regions on either side of the line absorption re- 
gion and thus are independent of changes in the spectral 
slope due to reddening. There is no viable continuum 
region to the red of the CO band head, so instead we 
extrapolate a linear fit to the continuum blueward of the 
feature and integrate the depth of the absorption under 
this extrapolated continuum. Although the slope of the 
fitted continuum will change with reddening, the inte- 
grated absorption equivalent width under the continuum 
should stay constant. 

Table |31 gives the observed spectroscopic indices for 
each star, and when available, colors and magnitudes. 
The photometry is from Frogel. Kuchinski & Ticd3 
l|1995D for Filler 1,'Kuchinski ct al. (1995) for Terzan 2, 
and Frogel & Sarajedini (1998, private communication) 
for Terzan 4. Measurement uncertainties for the spectro- 
scopic indices are typically 0.3, 0.4, and 0.6 A for Na, Ga, 
and CO, respectively, calculated from the errors associ- 
ated with fitting the continuum. However, an analysis 
of duplicate measurements taken over several years and 
with different instruments yields slightly higher values 
(see Table 3 of Paper I) . 

Column 2 of Table |21 lists the stars observed in 
each cluster. The identifi c ations for Filler 1 are from 
iFrogel. Kuchinski fc Tiedd lll995fl. and th e identifica- 
tions for Terzan 2 are from iKuchinski et al.l (11995). The 
numbering schemes for the remaining clusters are com- 
pletely arbitrary. We therefore provide finding charts for 
NGC 6256 in Figure □ NGC 6539 in Figure HP 1 in 
Figure 13 Palomar 6 in Figure 01 and Terzan 4 in Fig- 
ure O The last column of Table |21 lists the metallicity 
derived for each star as discussed in the next section. 

For clusters without calibrated near-IR photometry 
(either our own or that of others), we constructed in- 
strumental color-magnitude diagrams from new JK im- 
ages taken under non-photometric conditions, and chose 
the brightest stars near the cluster centers with colors 
and magnitudes consistent with membership. For clus- 
ters that physically lie within particularly crowded bulge 
fields, this may still not preclude the inclusion of non- 
members, since bulge giants will suffer from the same 
degree of reddening and will have colors and magnitudes 
similar to those of giants in the clusters, especially the 
metal rich clusters. Table01gives the cluster mean values 
for the indices and photometric parameters, along with 
the star-to-star intra cluster dispersions. The two entries 
for HP 1 will be explained later. 

Figure El illustrates, for each cluster except HP 1, the 
average of the individual stellar spectra. Here we have 
indicated the three features measured, the Na i doublet, 
the Ca I triplet, and the first band head of CO. For ref- 
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Table 1. Clusters observed & some previous metallicity determina- 
tions 



Cluster 


I 


b 


E{B - VY 




[m/H]''w'85 


[Fe/H]d^28g 


[Z/ZqYb9S 


[Fe/H]''^f03 


NGC 6256 


-12.2 


3.3 


1.03 




-1.56 




-1.01 


-0.70 


NGC 6539 


20.8 


6.8 


0.97 


-0.66 


-1.05 






-0.66 


HP 1 


-2.6 


2.1 


0.74 




-1.68 


-0.56 


-1.09 


-1.55 


Liller 1 


-5.2 


-0.2 


3.06 


-0.21 


-0.29 


+0.20 


+0.08 


+0.22 


Palomar 6 


2.1 


1.8 


1.46 


-0.74 


+0.22 




-0.09 


-1.09 


Terzan 2 


-3.7 


2.3 


1.57 


-0.47 


-0.54 


-0.25 


-0.26 


-0.40 


Terzan 4 


-4.0 


1.3 


2.35 


-0.21 


-0.29 


-0.94 


-0.61 


-1.60 


A[Fe/H]' 








0.37 ±0.62 


0.23 ±0.62 


0.64 ±0.08 


0.50 ± 0.28 


0.14 + 0.45 



"from the 2003 version of the 'Harris ('1996') catalog 

'^from integrated spectroscopy (Zinn & West 1984il) or NIR photometry IMalkarill982l) 
'^from dereddened subgiant colors iWebbink 1 9851) 

"^from integrated measurements of the Ca II triplet i ^mandrofl^^ Zinnll 1 98'M) 
•^from integrated measurements of the Ca II triplet i Bica et alj|1998|) 
compared to our final values given in Table m 



Table 2. Equivalent Width Measurement Intervals 



Feature 


Wavelengths (/im) 




Na I feature 


2.2040 - 2.2107 


* 


Na I continuum 


2.1910 - 2.1966 




Na I continuum 


2.2125 - 2.2170 




Ca I feature 


2.2577 - 2.2692 




Ca I continuum 


2.2450 - 2.2560 




Ca I continuum 


2.2700 - 2.2720 




^^CO(2,0) band 


2.2915 - 2.3025 




^2CO(2,0) continuum 


2.1900 - 2.2010 




12CO(2,0) continuum 


2.2110 - 2.2220 


• 


12CO(2,0) continuum 


2.2330 - 2.2600 


12CO(2,0) continuum 


2.2680 - 2.2800 


• 


12CO(2,0) continuum 


2.2860 - 2.2910 





N 
4 



162 



.51 



93 



Fig. 1. — The 2' x 2' / band finder chart for the stars observed 
in NGC 6256. 



erence, we have also marked a strong Fe i line and the 
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Fig. 2. — The 2' x2' K band finder chart for the stars observed 
in NGC 6539. 



second band head of CO. 

Figure [7| shows the individual stellar spectra for all 
stars observed in HP 1. Note that the first two of the 
HP 1 stars have significantly weaker absorption features 
than the other four stars. For comparison, we also show 
the spectrum of a low-latitude Galactic bulge star with 
similar colors and ind ices as the "strong" HP 1 stars from 
iRamfrez et all HHdS). 

3. DERIVATION OF CLUSTER METALLICITIES 

Paper I investigated techniques to use spectroscopic 
features measured in medium-resolution infrared spec- 
tra to measure globular cluster metallicities. The three 
strongest features in the K band spectrum are the Na i 
doublet, the Ca i triplet, and the first band head of CO. 

It is well known that each of these features has its own 
problems as a metallicity indicator. [Na/Fe] can show 
significant star-to-star variations within a single clus- 
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Table 3. Near-IR Spectral Indices, Colors, and Luminosities 



EW (A) 



Cluster 


btar 




(J - K)o 


IN a 


Oa 


cu 


n?^ /xji c 
[rcf ti\ 


NGC 6256 


51^ 






0.60 


0.71 


6.52 


-1.37 


NGC 6256 


93 






0.67 


0.85 


5.51 


-1.37 


NGC 6256 


95 






0.79 


0.37 


3.37 


-1.41 


NGC 6256 


162 






1.01 


0.40 


4.20 


-1.31 


NGC 6256 


B 






0.71 


0.46 


8.05 


-1.29 


NGC 6539 


1 






1.99 


2.84 


16.83 


-0.67 


NGC 6539 


2 






1.88 


2.35 


15.27 


-0.76 


NGC 6539 


3 






1.85 


2.71 


16.45 


-0.72 


NGC 6539 


4 






1.31 


1.55 


12.34 


-0.99 


HP 1 


46 






0.63 


0.93 


5.73 


-1.38 


HP 1 


57 






4.80 


4.30 


22.69 


-0.25 


HP 1 


191 






4.63 


4.01 


21.35 


-0.29 


HP 1 


435 






0.62 


1.21 


11.92 


-1.22 


HP 1 


A 






4.03 


2.98 


21.14 


-0.42 


HP 1 


B 






2.83 


3.22 


18.50 


-0.48 


Liller 1 


6 


-6.50 


1.17 


3.57 


4.63 


21.43 


-0.16 


Liller 1 


7 


-6.45 


1.13 


2.35 


5.34 


14.53 


-0.59 


Liller 1 


157 


-6.80 


0.93 


4.22 


5.62 


21.23 


-0.23 


Liller 1 


158 


-6.62 


1.13 


3.15 


4.32 


20.93 


-0.32 


Liller 1 


160 


-6.50 


1.27 


2.22 


2.04 


12.28 


-0.66 


Liller 1 


162 


-6.05 


1.16 


2.46 


3.27 


18.00 


-0.45 


Liller 1 


166 


-5.47 


0.95 


3.07 


3.16 


11.58 


-0.56 


Liller 1 


299 


-6.73 


1.38 


3.93 


5.34 


21.88 


0.07 


Palomar 6 


lb 






3.81 


4.71 


21.42 


-0.16 


Palomar 6 


2 






2.65 


2.94 


18.01 


-0.54 


Palomar 6 


3'' 






0.97 


-1.14 


12.64 


-1.03 


Palomar 6 


4 






1.11 


-1.94 


17.02 


-0.85 


Palomar 6 


5 






4.64 


3.04 


22.40 


-0.43 


Palomar 6 


11 






4.55 


4.75 


18.54 


-0.16 


Palomar 6 


14 






2.78 


1.77 


17.33 


-0.63 


Terzan 2 


1 


-5.82 


0.82 


1.37 


2.18 


12.45 


-1.07 


Terzan 2 


2 


-5.78 


0.95 


1.02 


2.39 


17.15 


-0.83 


Terzan 2 


3 


-5.73 


0.99 


4.86 


4.03 


22.39 


0.27 


Terzan 2 


4 


-5.72 


0.85 


1.93 


2.36 


18.16 


-0.74 


Terzan 2 


5 


-5.46 


0.97 


0.83 


2.85 


16.67 


-0.79 


Terzan 2 


7 


-5.36 


0.88 


1.32 


3.23 


16.58 


-0.76 


Terzan 2 


8 


-5.33 


0.87 


0.82 


1.83 


11.67 


-1.04 


Terzan 4 


185 


-6.44 


0.74 


0.97 


0.91 


6.28 


-1.49 


Terzan 4 


333 


-6.25 


0.84 


1.16 


-0.15 


7.67 


-1.33 


Terzan 4 


399 


-6.69 


0.73 


-0.07 


1.96 


-0.78 


-1.90 


Terzan 4 


424 


-6.11 


0.69 


-0.01 


2.02 


1.45 


-1.77 


Terzan 4 


453 


-6.61 


0.85 


0.01 


0.85 


8.70 


-1.57 


Terzan 4 


461 


-5.92 


0.70 


-0.33 


0.31 


3.18 


-1.80 


Terzan 4 


474 


-6.22 


0.78 


0.35 


0.89 


9.20 


-1.48 



''The average of measurements taken on two nights. 
''Non-members based on the radial velocities of Le e fc Carne 

'^Based on the QQ3 solution for NGC 6256, NGC 6539, HP 1, and Palomar 6, while Liller 1 and Terzan 4 are based on the LLA9 solution, 
and Terzan 2 is based on the QQA9 solution (see Table l5t . 



Table 4. Cluster, Mean Values and Dispersions for Observed Parame- 
ters 



Cluster 


^ stars 


EW(Na) 


EW(Ca) 


EW(CO) 


Mko 


{J 


-K)o 


mean 


a 


mean 


(7 


mean 


fj 


mean 


a 


mean 


a 


NGC 6256 


5 


0.76 


0.16 


0.56 


0.21 


5.53 


1.86 










NGC 6539 


4 


1.76 


0.30 


2.36 


0.58 


15.22 


2.03 










HP 1 weak 


2 


0.63 


0.01 


1.07 


0.20 


8.83 


4.38 










HP 1 strong 


4 


4.07 


0.89 


3.63 


0.63 


20.92 


1.75 










Liller 1 


8 


3.12 


0.75 


4.21 


1.28 


17.73 


4.33 


-6.39 


0.44 


1.14 


0.15 


Palomar 6 


5 


3.15 


1.48 


2.11 


2.50 


18.66 


2.17 










Terzan 2 


7 


1.74 


1.43 


2.70 


0.74 


16.44 


3.60 


-5.60 


0.21 


0.90 


0.07 


Terzan 4 


7 


0.30 


0.56 


0.97 


0.80 


5.10 


3.86 


-6.32 


0.27 


0.76 


0.06 



Note. — 



Pal 6 Excludes the 2 known non-members. 
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(b). 



• • ^435 

« 



f 191 • 
(d) i . . . 



57 



1 • • 



^46 



4- ' * 



• 424 



Fig. 3. — The 1.5' x 1.5' / band finder chart for the stars 
observed in HP 1. The letters in parentheses mark stars observed 
spectroscopically bv iMinniti. Cl995a.) (a.b.c.d.cen). 



Fig. 5. — The 1' x 1' / band finder chart for the stars observed 
in Terzan 4. 



(c) 



i 



r 

3(b) 



r 

1(f) 



E-*- 




2.20 



Fig. 4. — The 2' x2' K band finder chart for the stars observed 
in Palomar 6. The letters in par e nthese s mark stars observed 
spectro scopically bv ILee CarnevI 120021) (a-g), and bv .Minnitii 
11995al) ('b.c.d.f'l. 



Fig. 6. — Average of the individual stellar spectra for each 
cluster observed, except for HP 1 (see Fig. |7J. The positions 
of the Na I doublet, a strong Fe I line, the Ca I triplet, and the 
first two band heads of CO are indicated. Only the ^-^00(2,0) 
band is used for determining EW(CO). Note the relative absence 
of spectral features in NGC 6256 and Terzan 4. The high S/N of 
the spectra is also evident. 



ter; iSneden. Pilachowski. fc Krafti l)200Cl() have measured 
fT[Na/Fe]~ 0.25 dex in M15 and M92. Calcium is an 
a-element, and thus suffers from a-enhancement, which 
varies from cluster to clus ter a nd systematically with 
overall abundance (,Carnev..l99fit) . The first-overtone CO 
band heads in the K band saturate at high metallic- 
ities, making them less sensitive to carbon abundance 
and more dependent on the microturbulent velocity. 

While each of these features has its own shortcomings 
as a metallicity indicator, taken together and averaged 
over several stars in a cluster, we have demonstrated that 



they nonetheless provide a precise metallicity determina- 
tion (Paper I). Including stellar luminosities and colors, 
which help constrain the stellar temperature and gravity, 
improves the technique even further. 

Taking advantage of these facts, in Paper I we derived 
four equations with which one can calculate [Fe/H] for 
globular cluster giants. These equations use either lin- 
ear or quadratic combinations of the measured equivalent 
widths of Na, Ca, and CO, and, optionally, the infrared 
color and luminosity (Table 10 of Paper I). These equa- 
tions are applied to each observed star in a cluster to 
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Na I Fe I Ca I '=00(2,0) '=00(3,1) - 

II III I IT 




2.20 2.25 2.30 

Fig. 7. — K band spectra of the six stars observed in HP 1, and 
the spectrum of an inner galactic bulge star with similar colors 
and luminosity for comparison I Ramirez et al. 2000). The first 
two stars, 46 and 435, clearly have much weaker features than the 
other four stars in HP 1. These first two stars are therefore referred 
to as the "weak" group, while the latter four are referred to as the 
"strong" group. 



predict [Fe/H] for that star. The unweighted average of 
these values is then the [Fe/H] assigned to the cluster. 
Solutions LL3 and QQ3, linear and quadratic, respec- 
tively, use only the three spectroscopic indices; there- 
fore, they are independent of reddening and distance to 
the clusters. Solutions LLA9 and QQA9 also incorporate 
the (J — K)o colors and Mk magnitudes. Thus, while 
[Fe/H] estimates from LLA9 and QQA9 use more infor- 
mation and can be more precise, they are also affected 
by uncertainties in the reddening and distance estimates. 
Table lists the coefficients for the four equations; it is 
an abbreviated version of Table 10 from Paper I. The last 
line of Tabic [S] is our estimate of the typical uncertainty 
of the metallicity derived from each solution. It is based 
on the scatter of the individual stellar determinations in 
each of the calibrating clusters with six to nine stars ob- 
served per cluster. It does not include any uncertainties 
introduced by the observations of clusters to which the 
calibration is being applied, e.g. errors in calculating the 
distance modulus or reddening for the solutions that em- 
ploy colors and absolute magnitudes. Nor does it include 
the absolute uncertainty of the metallicity scale, which 
is discussed in section 

Table lists the average abundances calculated for 
each cluster from the four solutions in Table El The 
second column indicates the number of stars observed 
in each cluster. The sigma is the star-to-star scatter for 
each cluster. While there are no significant differences 
in the means between the [Fe/H] values from any of the 
solutions, based on the estimated error and the observed 
star to star scatter for each solution, we choose a "best" 



metallicity for each cluster and give the results in Ta- 
ble |HI These are computed on a cluster- by-cluster basis, 
as described in section 0] The metallicity of each star, 
as used in the calculation of the "best" value for each 
cluster, is also listed in the last column of Table 01 

3.1. How many Stars are Enough? 

One may ask how many stars are required to obtain 
a reliable measure of the metallicity of a cluster. In the 
clusters presented here we have tried to obtain spectra 
of more than 5 stars in each; however, because of time 
or weather constraints, this may not always have been 
possible; for example, in NGC 6539 we only have spec- 
tra of four stars. In such cases, where only a few stars 
have been measured, it becomes desirable to know with 
what accuracy and precision one can determine the clus- 
ter metallicity. Is it possible, for example, to get away 
with only measuring one or two or three stars with our 
technique? 

One important consideration is the potential contam- 
ination from superposed field stars. With a large sam- 
ple of spectra it is possible to effectively reject a few 
non-members if their metallicities are sufficiently differ- 
ent from that of the cluster. However, with a reduced 
sample size it becomes impossible to distinguish between 
members and nonmembers. 

A way to test the robustness of our metallicity determi- 
nations against small sample sizes is to run Monte Carlo 
simulations on some of the clusters that have larger sam- 
ples, and therefore relatively well-determined metallici- 
ties. To simulate the results of only measuring N stars 
in a cluster, we randomly select stars from the observed 
sample and average their metallicities to determine that 
of the cluster. By repeating this process 100, 000 times, 
we accumulate the distribution of errors, and can there- 
fore estimate the accuracy and precision with which one 
can determine the cluster metallicity using only a sample 
of N stars. 

Two illustrative examples are the spectra of the clus- 
ters Terzan 2 and Terzan 4. In Terzan 4 we measured 
seven stars that have a reasonably narrow distribution of 
metallicities; therefore, picking even a few stars gives a 
reasonable result for the cluster. On average, only mea- 
suring one star will give the correct answer, but the dis- 
tribution of results has a sigma of 0.19 dex. Averaging 
the results from two stars, the width of the distribution 
is reduced to 0.12 dex, and with three stars it drops to 
0.09 dex. 

Terzan 2 also has a relatively narrow distribution of 
metallicities, except that, of the seven measured stars, 
one appears to be a field star with a much higher metal- 
licity. In our analysis in Sec 14.61 we recognize this fact, 
and disregard the outlier when determining our "best" 
value for the cluster. However, with a smaller sample, it 
would be much more difficult to recognize contaminants. 
This is shown in the simulations, where on average, the 
deviation from the "best" value is 0.16 dex. Only mea- 
suring one star yields a distribution with a sigma of 0.42 
dex; the mean of two stars give a sigma of 0.27 dex, three 
a width of 0.20 dex, etc. 

Thus, the appropriate number of spectra for an ac- 
curate measure of the cluster metallicity (where accu- 
rate means to within the accuracy of the technique) ap- 
pears to be most dependant on the severity of the non- 
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LL3 LLA9 QQ3 QQA9 



variable 


coef 




coef 


(7 


coef 




coef 


cr 


, 

corife t 


-1.663 


0.057 


-1.451 


0.18 


-1.811 


0.074 


1.097 


1.16 


EW(Na) 


0.182 


0.029 


0.202 


0.036 


0.389 


0.065 


0.130 


0.091 


EW(Na)2 










-0.047 


0.013 


0.016 


0.020 


EW(Ca) 


0.057 


0.025 


0.025 


0.026 


-0.030 


0.051 


0.058 


0.056 


EW(Ca)2 










0.024 


0.012 


-0.004 


0.014 


EW(CO) 


0.027 


0.005 


0.026 


0.005 


0.043 


0.013 


0.028 


0.016 


EW(CO)2 










-0.001 


0.000 


0.000 


0.001 


{J-K)q 






0.749 


0.23 






5.240 


1.84 


{J-K)l 














-2.313 


0.93 


Mk 






0.151 


0.033 






1.812 


0.40 
















0.147 


0.034 


a (est) 


0.11 




0.07 




0.09 




0.06 





Table 6. Cluster Abundances 







LL3 




LLA9 


QQ3 


QQA9 


Cluster 


A^ stars 


[Fe/H] 


cr 


[Fe/H] 


a 


[Fe/HJ 


cr 


[Fe/H] 


a 


NGC 6256 


5 


-1.34 


0.04 






-1.35 


0.05 






NGC 6539 


4 


-0.80 


0.14 






-0.79 


0.14 






HP 1 weak 


2 


-1.25 


0.13 






-1.30 


0.11 






HP 1 strong 


4 


-0.15 


0.24 






-0.36 


0.11 






Liller 1 


8 


-0.38 


0.29 


-0.36 


0.25 


-0.31 


0.28 


-0.30 


0.29 


Palomar 6 


5 


-0.47 


0.44 






-0.52 


0.25 






Terzan 2 


7 


-0.75 


0.38 


-0.77 


0.40 


-0.83 


0.27 


-0.71 


0.45 


Terzan 4 


7 


-1.42 


0.15 


-1.62 


0.21 


-1.53 


0.28 


-1.60 


0.23 



Note. — Pal 6 Excludes the 2 known non-members. 



member contamination. If one expects no contaminants, 
e.g. an isolated cluster, three stars are sufficient to re- 
duce the stochastic noise inherent when using spectra of 
the signal-to-noise ratio (S/N) presented here (consistent 
with what was found in Paper I). However, in more pol- 
luted regions, such as the Galactic disk or bulge, more 
stars are required to be able to effectively reject non- 
members. Furthermore, the number of stars required in- 
creases as the difference between the cluster metallicity 
and the contaminating population decreases. 

3.2. The Metallicity Scale 

The calibration of the current technique used the 
met alhcities listed in the 1999 version of the Harris cat- 
alog l)Harrislll99fiD of Galactic globular clusters, which 
is a continually changing average of parameters taken 
from literature. Thus, while Harris started out on the 
[zinn & West ( 1984) scale, the catalog has been evolv- 
ing as more and more values from high-resolution spec- 
troscopy are incorporated. However there exist signif- 
icant discrepancies observed between metallicities ob- 
tained with various techniques, and on the absolute cal- 
ibratio n of the meta llicity scale itself. 

The iZinn fc Wesl l|1984 hereafter ZW84) metaUicity 
scale was one of the first such metallicity scales. Cal- 
ibrate! primarily from the spectroscopy of Cohen (see 
lCohenl[l983l and references therein), it used photomet- 
ric indices measured from cluster integrated light, and 
was applied to a total of 121 clusters. This scale has 
been updated bv lRutledee. Hesser fc StetsonI 1)19971) us- 
ing low-resolution spectroscopic measurements of the red 



Ca II t riplet in individual gi a nts. 

The ICarretta fc GrattonI llT997i hereafter CG97) 
[Fe/H] scale is based on a uniform analysis of high- 
resolution spectra of more than 160 bright red giants 
in 24 globular clusters. Figure 5 of CG97 shows a 
comparison between their system and that of ZW84: 
a nonlinear relation where the ZW84 values are ~ 0.1 
dex too high at the metal-rich end ([Fe/H]> —1), and 
~ 0.2 dex too low for clusters of intermediate metallicity 
(-2 <[Fe/H]< -1). 

The recent work by iKraft fc Ivan^ l|200.'l hereafter 
KI03), using measurements of Fe II in high-resolution 
spectra of giants in 16 clusters, confirms that the ZW84 
[Fe/H] scale is nonlinear with respect to the true [Fe/H] 
scale. Their Figure 3 shows the clearly nonlinear rela- 
tion between theirs and the ZW84 scale, regardless of 
the model atmosphere they choose. However, Figure 2 
of KI03 illustrates the comparison of their scale with 
that of CG97, showing that the CG97 values are too 
high, with an offset of ^ 0.2 dex at the lowest metal- 
licities ([Fe/H]= —2.4), and agreement at the highest 
([Fe/H]= -0.7). 

Since our calibration is not strictly based on any one 
of these three different metallicity scales, we have per- 
formed a comparison of the metallicities of the 15 clusters 
used to calibrate our technique. The results of this com- 
parison are shown in Figure |S1 This shows that at low 
metallicities, CG97 gives values that tend to be higher 
than our calibration, and KI03 gives values that tend to 
be lower than our calibration. At higher metallicities, we 
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Table 7. Comparison of Calibrating Clusters 



Scale 


N 


A [Fc/H] 


cr 


ZW84 


15 


+0.013 


0.123 


Rut97 


12 


-0.014 


0.122 


CG97 


6 


+0.090 


0.121 


KI03 


9 


-0.076 


0.110 



Table 8. Best Values 



Cluster 


[Fe/H] 


NGC 6256 


-1.35 ± 0.20 


NGC 6539 


-0.79 ± 0.20 


HP 1 


-1.30 + 0.20 


Liller 1 


-0.36 + 0.20 


Palomar 6 


-0.52 + 0.20 


Terzan 2 


-0.87 + 0.20 


Terzan 4 


-1.62 + 0.20 



are lower than ZW84, but not as low as KI03. 
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Fig. 8. — Comparison of our calibration globular cluster 
metallicities as derived on different scales, as a function of the 
metallicity used in the calibration of our infrared spectroscopic 
technique. The difference is in the sense that if a cluster has 
a higher metallicity compared with what we used in our cali- 
bration, it will show a positive A [Fe/H 1. The crosses represent 
the clusters listed in Zinn & West" (1984) , the triangles are clus- 
ters from Rutlcdgc, Hesscr & Stetson (1997), the c ircles are from 
ICarretta &: GrattonI 119971) . and the stars are from lKraft fc Ivand 
120031) . 

Table [7| summarizes the data in Figure |S1 it includes 
the metallicity scale, the number of calibrating clusters 
measured on each system, the mean difference in metal- 
licity, and the dispersion of the difference in columns (1) 
through (4) respectively. In the mean, our values lie 
closest to those of ZW84 and Rut97, while they are in- 
termediate between those of CG97 and KI03. 

Thus, although we can reproduce the metallicities of 
our calibration stars to better then ±0.1 dex, in Table|Sl 
we formally quote our uncertainty as ±0.2 dex to remind 
the reader that this technique is only as good as the 
metallicity scale of the calibrating clusters. 

4. NOTES ON INDIVIDUAL CLUSTERS 
4.1. NGC 6256 



iBica et all l|1998j) have determined [Z/Zq] for a num- 
ber of metal rich clusters based on measurements of the 
infrared Ca ii triplet in integrated spectra. However, due 
to a-element enhancement in metal-rich clusters, their Z 
values are expected to be greater than the true [Fe/H]. 
For NGC 6256 they derive [Z/Zq]= -1.01 ± 0.15, but 
also note that this may only be an upper limit due to 
contamination of the integrated spectrum by metal-rich 
b ulge stars. 

iQrtolani. Barbuv fc Bical l)1999tl measured the optical 
CMD of NGC 6256, and noted that it appears to have 
an extended blue horizontal branch and an RGB similar 
to that of NGC 6752 ([Fe/H]= -1.54). Based on these 
features, and taking into account the high Bica et al. 
value, they estimate [Fe/Hjw —1.3. 

Using infrared spectra of five giants, we apply the 
metallicity solutions based on only the equivalent widths 
as there is no calibrated photometry for this cluster. LL3 
gives —1.34 with a dispersion of 0.04, and QQ3 gives 
— 1.35 with a dispersion of 0.05. In this case the er- 
ror of the mean metallicity will not be set by the dis- 
persion of observed values (a/7V^/^), but rather by the 
accuracy of the metallicity solution (see Table [SJ. We 
therefore choose the QQ3 solution, which has a slightly 
smaller estimated error, and assign [Fe/H]= — 1.35±0.09 
to NGC 6256. 

Ortolani's estimate of [Fe/H]w —1.3 is in excellent 
agreement with our value, and the Bica et al. measure- 
ment is also within range of our value given their quoted 
caveats (i.e. that their value of —1.01 is an upper limit). 
However, the value of —0.70 listed in the Harris catalog 
appears to be in need of revision downward. 

4.2. NGC 6539 

The only two metallicity determinations for NGC 6539 
are those of ^nn & West (1984), who found [Fe/H] = 
—0.66 ± -15 ba sed on integrated optical spectra, and 
iWebbinkI l(T985|l . who estimated [m/H]= -1.05 ± 0.28 
from dereddened subgiant colors. 

Our new measurement is based on spectra of four stars 
in NGC 6539. As there is no calibrated infrared photom- 
etry, we measure metallicities using only the equivalent 
widths. The LL3 solution yields —0.80 with a dispersion 
of 0.14, and the QQ3 solution gives —0.79 with a dis- 
persion of 0.14. These star-to-star dispersions are very 
small, and the error of the mean metallicity is therefore 
set by the estimated accuracy of the technique. Tak- 
ing the QQ3 value, we assign [Fe/H]= -0.79 ± 0.09 to 
NGC 6539. 

Both the Zinn & West and Webbink measurements 
agree with our value, as well as the Harris catalog value 
of -0.66. 

4.3. HP 1 

IWebbinkI l|1985H first estimated the metallicity of HP 1 
to be [m/H]= —1.68 ± 0.28 based on dereddened sub- 
giant colors. Soon thereafter lArmandroff fc ZinnI l)1988f) 

found the metallicity to be more than 10 times higher, 
[Fe/H]= —0.56 ± 0.12, from the strength of the infrared 
Ca II triplet in an integrated spectrum of HP 1. 'MinnitiJ 
(^95a) measured an even higher value of [Fe/H] = 
—0.30 ± 0.20 based on the strength of several Mg and 
Fe features in six bright cluster giants (indicated in Fig- 
ure chosen on the basis of their position in the IR 
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CMD. Unfortunately, there is no overlap between Min- 
niti's sample and our own; however, it is worthwhile to 
note that based on his radial velocities, stars c and d are 
non-members, while a and b have velocities similar to 
t hat of the central grouping of star s ("cen"). 

lOrtolani. Bica fc BarbuvTljlQQTal) claim that these very 
high metallicities are due to contamination by metal-rich 
bulge giants. Using VI CMDs they find that HP 1 has 
an extended blue HB and a CMD that resembles that of 
NGC 6752, and they estimate [Fe/H]« -1.5 for HP 1. 

iBica et al. ( 1998, see note under NGC 6256) agree that 
contamination from metal-rich bulge giants is very im- 
portant in HP 1. Their integrated spectrum of the cen- 
tral 10" yields a low metalHcity of [Z/Zq]^ -1.09±0.15, 
while the spectrum covering the whole cluster yields a 
much higher metallicity of [Z/Zq]~ —0.37. 

Usin g a CO filter to isolate cluster members, iDavidgd 
Ij2000l) found that the RGB of HP 1 is very similar to that 
of M 13 in a if, J - K CMD, implying [Fe/H] 1.6. 

Our infrared spectroscopic analysis uses six stars mea- 
sured in HP 1. As mentioned earlier, we have divided 
the HP 1 spectra into "weak" and "strong" groups based 
on the strengths of their absorption features (see Fig- 
ure O . Identification of these stars in the Ortolani et 
al. CMD shows that stars 57, 191, and A are very red, 
with V — I > 4, and lie on the curved RGB of the metal- 
rich bulge. Stars 46, 435, and B are not as red, with 
V" — / < 3, and appear to lie on the more vertical RGB 
of HP 1. Therefore, we conclude that the two "weak" 
stars, 46 and 435, are the true HP 1 members, while the 
other "strong" stars are surrounding bulge stars. 

Based on the weak group, consisting of only stars 46 
and 435, we estimate the cluster metallicity using only 
the measured equivalent widths: —1.25 with a dispersion 
of 0.13 from the LL3 solution, and —1.30 with a disper- 
sion of 0.11 from the QQ3 solution. The dispersions are 
so small that the error of the mean is determined by the 
accuracy of the technique, and we assign a metallicity of 
[Fe/H]= -1.30 ± 0.09 for HP 1. 

It seems clear that HP 1 has a low metallicity, and 
higher estimates were due to contamination by metal- 
rich bulge stars. The measurements of Ortolani et al.(« 

-1.5), Bica et al. (-1.09 ± 0.15), and Davidge { 1.6), 

as well as the Harris catalog value of —1.55, are all in 
reasonable agreement with our new measurement. 

4.4. Liller 1 

iMalkanI l)1982D first estimated the metallicity of Liller 1 
to be just below solar, [Fe/H]= —0.21, b ased on inte- 
grated near-IR photometry (also listed in "Zi nn fc Westi 
11284; Zinn 1985). Armandroff & Zinn (1988) later raised 
the value to above solar, [Fe/H]= -t-0.20±0.12, from their 
measurement of the integrated Ca ii triplet. 

Studies based on the CMD morphology seemed to 
confirm the near-solar, or even supersolar metallic- 
ity of Liller 1. The optic al RGB measured by 
lOrtolani. Bica fc Barbuvl l)1996f) shows strong curvature 
and indicates that Liller 1 is considerably more metal- 
rich than NGC 6528, i.e. [Fe/H]> -0.17. The slope of 
the giant branch in the infrared also gives a high metal- 
licity, [Fe/H]= -hO.25 ± 0.30 (Froecl. Kuchinski & Ticdc 
11295'), although this is an extrapolation of an empirical 
relation that has only been calibrated for [Fe/H]< —0.3. 

The first suggestion of a significantly sub-solar metal- 



licity came from medium-resolution H and K band spec- 
tra of the central 4 .4" x 6.6" of Liller 1 obtained by 
lOridia et all l)1997l) . They calculated [Fe/H] for the 
integrated light based on the strength of the CO (6,3) 
^ absorption feature at 1.62 /im an d found [Fe/H] = 
-0.29 ± 0.3. Although soon thereafter iBica et all l|199a[ 
see note under NGC 6256) derived [Z/Zq\= -h0.08±0.15 
for this cluster based on the integrated Ca li triplet. 

Most recently, Davidge (2000) has estimated that the 
metallicity of Liller 1 is comparable to, or slightly less 
than that of NGC 6528 ([Fe/H]< -0.17) based on filter 
determined CO band s t rength s and infrared colors, and 
lOriglia. Rich fc Castrol l)2002D have used high-resolution 
infrared H band spectra of two bright giants to measure 
an abundance of [Fe/H]= -0.3 ± 0.2. 

We have obtained spectra of eight stars in Liller 1. 
Using only the equivalent widths, we find a metallic- 
ity of —0.38 with a dispersion of 0.29 with LL3 and 
—0.31 with a dispersion of 0.28 with QQ3. How- 
ever, including the stellar luminosities and colors from 
Frogel, Kuchinski fc Tiedc ( 199 5.), assuming (m — Af )n = 
14.91 and E{B - V) = 3.06 ijHarrid Il996|) . we obtain 
—0.36 with a dispersion of 0.25 with LLA9 and —0.30 
with a dispersion of 0.29 with QQA9. These values 
are all consistent with one another, and we choose the 
LLA9 solution, which has the smallest dispersion, giving 
[Fc/H]= —0.36 ± 0.09 as our best metallicity measure- 
ment for Liller 1. 

The recent medium- and high-resolution spectroscopic 
measurements of individual giants by Origlia et al. are 
in perfect agreement with our value, although the esti- 
mates based on the CMD morphology or integrated spec- 
tra seem to give values which are too high. Harris's value 
of +0.22 also needs revision downwards. 

4.5. Palomar 6 

The first estimate of the metallicity of Pa lomar 6 
came from the integrated infrared photometry of lMalkai] 
J 19821), who found [Fe/H] = -0- 74. Over a decade later, 
Ortolani, Bica fc Barbuvf (|1995|) published the first opti- 
cal CMD and estimated that its metallicity is compara- 
ble to NGC 6356, i.e. [Fe/H] 0.4, based on the GB 

morphology. 

The highest metallicity estimates come from iMinnitil 
()1995a|) and ^ica et al. (1998). Minniti measured 
[Fe/H] =-1-0.20 ± 0.30 based on optical spectroscopy of 
four giants selected via IR photometry and marked with 
letters b,c,d,f in Figure^ Bica et al. found [Z/Zq\ = 
—0.09 ± 0.15 based on integrated spectra of the infrared 
Ca II triplet. 

Most recently, iLee fc CarnevI l|2002f) measured a very 
low [Fe/H]= —1.19 ± 0.18 based on the slope of the in- 
frared GB, and — 1.08±0.06 from high-resolution IR spec- 
tra of three giants. They also measured radial velocities 
from high-resolution IR echelle spectra of seven giants, 
marked with letters a-f in Figure 01 Their radial veloci- 
ties indicate that star 2(g) is a member, while stars 1(f) 
and 3(b) are nonmembers. We thus exclude stars 1 and 
3 from our analysis of the cluster metallicity. 

^ They point out that the CO bands in the H band are better 
indicators of [Fe/H] because they are not as heavily saturated as 
those in the K band, and are thus more sensitive to abundance 
and less sensitive to the microturbulent velocity. 
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As no calibrated infrared photometry is available for 
Pal 6, our metallicity determinations are based on the so- 
lutions that use only the equivalent width measurements. 
With the five stars measured in Palomar 6, we find —0.47 
with a dispersion of 0.44 using LL3 and —0.52 with a 
dispersion of 0.25 using QQ3. We choose the quadratic 
solution for our best value, as it has the smallest disper- 
sion among the individual stars, assigning a metallicity 
of [Fc/H]= -0.52 ± 0.11 to Palomar 6. 

This cluster has an exceptionally large range of metal- 
licity measurements, ranging from —1.2 to -1-0.2, where 
the Harris cata log has recently adopt ed a value of —1.09, 
primarily from lLee fc Carnevl l)2002fl . Our measurement 
of —0.52 ± 0.11 falls in the middle of t his range, but 
is much higher than that measured by iLee fc Carnevl 
l|2'002'l . One possibility is that more of our stars are 
non-members, belonging instead to the metal-rich bulge 
population. If we make the assumption that stars 5 and 
11, which have the strongest absorption features (see Ta- 
bleOl, s-lso non-members, keeping only stars 2, 4, and 
14, the resulting metallicity would be —0 .7 ± 0.2 (LL3), 
nearly consistent with that obtained by iLee fc Carnevl 
lf20f)l . 

4.6. Terzan 2 

iMalkanI l)1982|) first estimated the metallicity of 
Terzan 2 to be [Fe/H]= —0.47 from integr ated infrared 
photometry (see also iZinn fc Wesl 119841: iZinnI Il985j) . 

The ensuing two studies both foun d Terzan 2 to have a 
slightly higher metallicity of —0.25. lArmandroff fc ZinnI 
(jiaksj measured [Fe/H]= -0.25 ± 0.12 from integrated 
spectroscopy of the Ca ii triplet. iKuchinski et al.l l|l995D 
measured [Fe/H]= — 0.25±0.25 based on the slope of the 
R GB in their infrared CMP 

lOrtolani. Bica fc Barbuvl l)1997 bf) estimated a slightly 
lower metallicity from their optical CMD. Measuring 
the difference in V magnitude between the HB and the 
brightest part of the RGB, they concluded that the 
metallicity is between that of 47 Tuc and NGC 6356, 
or [Fe/H] 0.55. 

The most recent measurement is that of iBica et ahl 
(jl998l) . who find [ZIZq]= -0.26 ± 0.15 from integrated 
optical spectra (see note under NGC 6256). 

From our infrared spectra, using only the equivalent 
widths measured in seven stars, we find a metallicity of 
—0.75 with a dispersion of 0.38 (LL3) and —0.83 with 
a dispersion of CK27 (QQ3). Inc l uding the luminosities 
and colors from IKuchinski et al.l l|1995l), a ssuming (m — 
M)o = 14.70 and E{B - V) = 2.35 (garris 19961, we 
find -0.77 with a dispersion of 0.40 (LLA9) and -0.71 
with a dispersion of 0.45 (QQA9). 

Looking at the distribution of individual stellar metal- 
licities (last column of Table I^J, star 3 is a > 2fT outlier, 
regardless of which solution is employed. If we disregard 
star 3, assuming that it is a non-member, the dispersion 
of the metallicities is reduced significantly, and the error 
of each solution becomes the limiting factor determin- 
ing the uncertainty in the metallicity of the cluster. We 
therefore choose the QQA9 solution as the "best" and 
assign a metallicity of [Fe/II]= -0.87 ±0.06 to Terzan 2. 

Our measurement is the lowest to date, significantly 
lower than the Harris catalog value of —0.40. Only the 
work of Ortolani et al. is consistent with our determina- 
tion. 



4.7. Terzan 4 

Because of its crowded, metal-rich environment of 
the Galactic bulge, integrated studies of Terzan 4 tend 
to yield artificially high met allicities . Ba. sed on in- 
tegrated infrared photometry IMalkanI l)1982f) measured 
[Fe/H]= —0.21. Using integrated spectroscopy of the in- 
frared Ca II triplet Armandroff & Zinn ( 1988) measured 
[Fe/H]= —0.94 ± 0. 12, a nd from the integrated spectra 
of Ter 4 IBica et al.l (IT998.^ found [Z/Zq]= -0.61 ± 0.15. 

However, studying th e opti cal CMD of Terzan 4, 
lOrtolani. Barbuv fc Bical l)1997D point out its similarity 
to M30, in particular the presence of a blue horizontal 
branch, which implies that Terzan 4 could be as metal- 
poor as [Fe/H]~ —2. 

Our spectroscopic analysis, using only the equiva- 
lent widths of the seven measured stars, yields metal- 
licities of —1.42 with a dispersion of 0.15 (LL3) and 
— 1.53 with a dispersion of 0.28 (QQ3). Including col- 
ors and magnitudes from Frogel fc Sarajedini (1998, pri- 
vate communication), assuming a distance modulus of 
(m-M)n = 1 4.795 and an extinction ofE{B-V) = 2.35 
l)Harriairi996D . we find —1.62 with a dispersion of 0.21 
(LLA9) and -1.60 with a dispersion of 0.23 (QQA9). 
The error of the mean of the LL3 solution is limited by 
the accuracy of the solution (0.11 dex), and therefore 
the LLA9 value gives the smallest error, and we assign 
[Fe/H]= -1.62 ± 0.08 to Terzan 4. 

The Harris catalog value of —1.60 is in perfect agree- 
ment with our new measurement, and the Ortolani et al. 
estimate of > —2 is also consistent. 

5. Lj CENTAURI 

It has been long known that u Centauri is unique 
among Galactic globular clusters in that it exhibits 
a substantial spread in stellar abundances. This 
was first seen in the breadth of the giant branch in 
early optical llWoollevI [l966: Dickens & WooUcv 1961 
iCannon fc Stohielll973t ILlovd Evans.. 1977.1 and infrared 
ijPersson et all 1198(H) color-magnitude diagrams. Later 
spectroscopic studies of individual giants fCohen"19811 
jMaUia fc Pagcl 1981; Gratton 1982; Pahoglou & Norril 
1989; Brown & Wallcrstcin Il993t iNorris fc Da Co^ 
il995; .Suntzeff fc Kraft,.1996^) showed the range to be ap- 
proximately one dex. The most recent photometric and 
spectroscopic work shows evidence for three subpopula- 
tions in uj Gen. The largest is a metal-poor group with 
a metallicity peak at [Fe/H] ^ —1.6, a moderately pop- 
ulated intermediate-metalhcity group at [Fe/H]~ —1.2, 
and a small number of metal-rich stars with [Fe/H]< 
—0.5, w hich make up about ^ 5% of the population 
llNorris. Freem an, fc Mighell 1996; Pancino et al. 200^ 
lOriglia et a,l]i2nol . 

We have obtained infrared spectra of 12 uj Centauri 
iants , selected from the photometry of iPersson et alJ 
19801) . Figure El shows the individual spectra, indicat- 
ing the wavelengths of some of the strongest absorption 
features. A mere visual inspection of these spectra re- 
veals that most are quite metal-rich, while V162 is quite 
metal-poor. 

Ta ble 01 li sts the infrared photometry from 
Perss oiTet all l)198nj) a ssuming^ (m — M)o — 13.62 
and E{B - V) = 0.12 ijHarris 19961 measured equiv- 
alent widths, and derived metallicities using each of 
the four equations from Paper I. We note that several 
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Fig, 9. — Spectra of the 12 lu Centauri giants. V164 is the 
average of two spectra taken on different nights. 



of our stars appear to fall into the "anomalously 
metal-rich" population, which only accounts for ^ 5% of 
cluster |Origfia_et al. 2003), while V162 is the only star 
measured from the predominant metal-poor population. 

Previous metallicity determinations of the stars mea- 
sured here a re also listed in the last t hree columns of 
Table El The lNorris fc Da Costal l(T995|) values were de- 
ter mined from high resoluti on optical echelle spectra. 
The ISuntzeff fc Krai llT99l values are based on mea- 
surements of the red Ca ii triplet calibrated from the 
line strengths measured in Galactic globulars (ZW cali- 
bration; an alternate calibration is based on 24 stars in 
common with Norris & Da Costa (1995)). The "other" 
values are explained in the table notes. 

On average, there is reasonable agreement between 
our values and those previously published. Taking 
the average of the four IR solutions listed in Table EJ 
the mean differenc e for the six stars in common with 
INorris k D a Costa' ( 1995) is -0.15 ± 0.27 dex. Compar- 
ing with ISuntzeff fc Kraft (,1996,1 . the mean difference is 
either 0.02 ±0.19 or —0.16 ±0.11, depending on whether 
we use their ZW or NDC calibration, respectively. 

However, when comparing with the four "other" metal- 
licity determinations listed in the last column of Table El 
most of which are from high-resolution spectroscopy, we 
seem to measure systematically higher values; the mean 
difference is 0.40 ± 0.13 dex. At this time it is unclear 
what could be causing this offset, although it is most 
likely not due to a-enhancement, since the stars in lu 
Cen have normal g-e nhancement compar ed with other 
globular clusters fe.g. lPancino 'erani200l . 



Looking at the observed relations of the measured 
equivalent widths in uj Centauri, for example EW(CO) 
versus EW(Na), EW(CO) versus EW(Na), and EW(Ca) 
versus EW(Na), we find that they are completely con- 
sistent with what was seen in the calibrating clusters 
of Paper I (see for example Figure 6 in Paper I). This 
indicates that, at least to within our measurement ac- 
curacy, the abundance ratios seen in these co Cen gi- 
ants are not vastly different from what we measured 
for our calibrating Galactic clusters. This is consistent 
with what has been measured for v arious elements using 
high-resolution sp ectroscopy (e.g. iPancino et al.l 120021: 
lOriglia et a,lJl200l . 

5.1. CO 

In Paper I we observed a tight relationship between the 
spectroscopic CO equivalent width and the photometric 
CO index. Taking the measurements of the photometric 
CO index from Pcrsson et al. (1980), wc find the same 
relation also holds true in to Centauri. 

Figure [TUl shows the spectroscopic EW(CO) as a func- 
tion of the photometric CO index for all the stars mea- 
sured in Paper I (dots) and the stars measured in tu Cen 
(open circles) . Both sets of data clearly follow the same 
relationship and arc well modeled by a simple linear least- 
squares fit. The solid line is the fit to only the data from 
Paper I, while the dashed line is to both data sets, except 
for the stars with overplotted x 's, which are greater than 
2.5(7 outliers and are not included in the fits. The fit that 
includes both data sets has slightly smaller errors than 
what was derived in Paper I, and is given in Equation^ 

EW{CO) = 2.97(±0.47) + 84.7(±3.4) xCO(mdea;) (1) 

6. DISCUSSION & CONCLUSIONS 

In this paper we present new metallicities for seven 
heavily reddened bulge globular clusters using medium- 
resolution infrared spectroscopy. Table El summarizes 
these results. These clusters are some of the most dif- 
ficult to study in the Galaxy, and comparing our new 
measurements with previous values, we find several sys- 
tematic trends as well as a few significant discrepancies. 

Figureim illustrates these differences and trends. Here 
we plot our new metallicity determinations from Table El 
with those from previous works listed in Table ^ The 
dots represent our measurements, where the uncertain- 
ties are taken to be ±0.2 dex, as discussed in Section rO 

The ijBica et al.lll998j) measurements are based on the 
strength of the red Ca ii triplet in the integrated cluster 
light. Ho wever, s ince Galactic globulars are a-enhanced 
(see lCarnevll996|, for a review of g-e nhancement in glob- 
ular clusters), the [Z/.Zq] values of iBica et all ifiggSL 
which measure an overall metallicity, will be slightly 
higher than the iron abundance, [Fe/H], Indeed, compar- 
ing with_ouTfinal_^undances, we find that in the mean, 
the lBica et al.l 1119981) values are systematically higher by 
0.50 ± 0.28 dex (see last row of T abled. 

The lArmandroff fc ZinnI l|1988L AZ88) measurements, 
also based on the strength of the red Ca il triplet, also 
give a large deviation from our final values: 0.64 ± 0.08 
dex. There are several possible explanations for this 
large discrepancy. AZ88 calibrated their summed Ca ii 
equivalent width (SVF) using a linear fit between TiW 
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Table 9. Observations and Analysis of u: Centauri Giants 
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84 


-5.78 


0.85 


1.14 


1.27 


17.56 


-0.91 


-0.97 


-0.98 


-0.96 


-1.36 


-1.16 


-1.34<= 


179 


-5.66 


0.92 


1.50 


1.64 


19.12 


-0.78 


-0.78 


-0.86 


-0.74 


-1.10 


-0.80 




201 


-5.71 


0.93 


1.25 


1.94 


21.80 


-0.74 


-0.75 


-0.90 


-0.69 


-0.85 






219 


-5.20 


0.85 


0.83 


1.21 


15.95 


-1.01 


-0.99 


-1.09 


-0.94 


-1.25 


-1.01 


-1.39*' 


248 


-5.42 


0.94 


0.89 


0.93 


9.39 


-1.19 


-1.12 


-1.19 


-1.08 


-0.78 


-1.03 




270 


-4.81 


0.79 


0.81 


1.33 


15.35 


-1.03 


-0.99 


-1.10 


-0.91 


-1.22 


-1.13 


-1.58S 


425 


-5.24 


0.98 


2.21 


2.45 


16.89 


-0.67 


-0.56 


-0.67 


-0.49 








447 


-5.13 


0.99 


2.53 


2.12 


14.56 


-0.69 


-0.54 


-0.67 


-0.47 




-0.25 




451 


-4.50 


0.87 


1.58 


2.01 


12.99 


-0.91 


-0.77 


-0.89 


-0.57 








513 


-5.08 


1.01 


2.02 


2.67 


14.51 


-0.75 


-0.61 


-0.71 


-0.52 






-0.90'^ 


V162 


-5.54 


0.75 


0.40 


0.81 


0.55 


-1.53 


-1.61 


-1.65 


-1.69 








V164' 


-5.44 


0.93 


1.61 


1.58 


12.80 


-0.93 


-0.88 


-0.91 


-0.85 









^Star IDs are from t he ROA catalog fWooiievllT96^ 
''Photometry is from 'Per sson et al.l 11980) 

' Norris & Da Costa 1 1995") from high resolution optical spectra. 

' Sunteeff & Kraft 1 1996), based on measurements of the Ca II triplet lines (ZW calibration). 
' Browi^^W allcrstcin 1 1993) from high resolution optical spectra. 
GrattonI ri982) from high resolution optical spectra. 
^[Z / ZQ]{ron\ Francois, Soite, & Spite 1 1988) from high resolution optical spectra. 
* 10riglia ct al.. 12003) from medium resolution infrared spectra. 
'V164 is the average of two spectra taken on different nights. 




0.1 

CO index (filters) 

Fig. 10. — Relationship between the spectroscopically deter- 
mined EW(CO) and the photometric CO index for the stars mea- 
sured in Paper I (dots) and uj Cen (open circles). The solid line 
shows the best-fit least-squares relation for the stars measured in 
Paper I only, and the dashed line is the best-fit relation including 
the stars measured in uj Cen, where in both cases stars with over- 
plotted x's are > 2.5(t outliers, which are not included in the fits 
(three of which are known large-amplitude variables). 
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seven clusters. However, it is now 
strength of the Ca ii hnes shows a non-hnear relation- 
ship with fFe/HlzivK4, especially for fFe/Hlzws4 > -0-7 
ijRutledge. Hesser fc StetsonI Il997t iKraft fc Ivaii3 120031 
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Fig. 11. — Comparison between previous metallicity determi- 
nations and our final [Fe/H] values determined in Section |2 The 
dots with error bars show our measurements (SF) and their ab- 
solute uncertainties. The previous measurements have been taken 
from Table although for simplicity we have omitted their er- 
ror bars. The vertical dashed line and shaded region illustrates 
the mean metallicity and the measured dispersion of the Galactic 
bulge (-0.21 ± 0.30: . Ramirez et al...200aV 



see also Section f3.2f) . This poor calibration could cer- 
tainly lead to overestimates of the metallicity; however, 
the most likely explanation is simply contamination from 
surrounding metal-rich bu lge star s . 

The Harris catalog ijHarrisI I1996|) IS an aver- 
age of globular cluster parameters taken from 
the literature. It is maintained on the web at 
http: / /physun. physics. mcmaster.ca/Globul ar. html. 
While originally on the lZinn & WestI l|1984l) scale, which 
appears slightly nonlinear and tends to overestimate the 
metallicities of the most metal-rich clusters, more recent 
versions of the catalog have incorporated many values 
fro m high-resolution spectrosc opy and are thus closer to 
the lCarretta fc GrattonI l)1997j) scale (see Section for 
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a discussion of the various metallicity scales, and how 
our technique compares). 

In general, the comparison of our values with those in 
the 2003 version of the Harris catalog is very close; the 
mean difference between the two sets of measurements is 
+0.14 ± 0.45 dex (last hne of Table [U, where the Har- 
ris values are slightly higher than ours. The explanation 
behind the higher metallicities found in previous studies 
seems to be primarily contamination by metal-rich bulge 
stars (see Figure ITT)) . The combination of very crowded 
fields, low-concentration clusters, and similar stellar lu- 
minosities makes it very difficult to distinguish between 
cluster members and superposed bulge stars. 

The largest discrepancies between our new values and 
those in the Harris catalog are those of NGC 6256 and 
Liller 1, which are listed with metallicities 0.65 and 0.59 
dex higher than we measure, respectively. Our lower 
values for both clusters are in good agreement with other 
recent studies. The Harris catalog also lists Palomar 6 
with a metallicity of —1.09, taken from the IR spectra 
of .Lee fc Carney (^2002), which is 0.57 dex lower than we 
measure. See Section 1^31 for a discussion of this cluster. 

As a result of this study we have reduced by three the 
number of Galactic globular clusters that have [Fe/H]> 
—0.7 according to the Harris catalog. The metallicity 
of NGC 6256 was lowered to -1.35 and NGC 6539 was 
decreased down to —0.79. Terzan 2, which lies within 10° 
of the Galactic center, was reduced to [Fe/H]=— 0.87. 

We have also observed 12 giants in the globular cluster 



Lo Centauri, all except one of which appear to lie in the 
intermediate- or metal-rich regions of oj Gen's tri-modal 
metallicity distribution. We compare our results for the 
eight stars that have previous metallicity determinations 
and find general agreement, although some of the most 
recent optical spectroscopy finds slightly lower metallic- 
ities than we predict. 
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